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; 1 INTRODUCTION 
•rH . 

, The puzzle of the origin of the X-ray-background (XRB) 
■ remains unsolved after over 30 years of study. Thermal 
C3 ' bremsstrahlung with a temperature of about 40 keV fits 
its spectrum between 3 and 60 keV very well (Marshall et 
al. 1980). Despite this agreement, observations of the cosmic 
microwave background (CMB) exclude the possibility that 
hot intergalactic gas is a major contributor to the XRB. 
The excellent fit of the CMB to the blackbody function 
strictly limits the amount of hot diffuse gas between the 
current epoch and the last scattering surface at redshift 
2 ~ 1000; such gas would distort the XRB spectrum through 
the Compton scattering. The upper limit to any such dis- 
tortion gives an upper limit to the contribution of the hot 
intergalactic gas to the XRB as < 3 per cent (Mather et 
al. 1990). 

It is therefore likely that the XRB is due to a contri- 
bution of different types of discrete sources. Active galactic 
nuclei (AGN) are likely candidates (see Fabian and Barcons 
1992 and references therein), especially since Seyferts and 
Quasars (QSO) provide a large fraction of the soft XRB be- 
low ~ 2 keV (Shanks et al. 1991; Hasinger et al. 1993). One 
of the major problems with an AGN origin of the XRB has 
been the apparent discrepancy of the XRB spectrum with 
the spectrum of resolved AGN. In the 2-10 keV band AGN 
spectra are too steep and subtraction of their contribution 
worsens the situation (this is the spectral paradox of Boldt 
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We discuss the possibility that a significant contribution of the hard X-ray Back- 
ground is the integrated emission from a population of galaxies undergoing advection- 
dominated accretion in their nuclei. Owing to poor coupling between ions and elec- 
trons and to efficient radiative cooling of the electrons, the accreting plasma is two- 
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there is some evolution associated with the spectral emissivity which must peak at a 
redshift ~ 2. We estimate that such galaxies contribute only to a small fraction of 
the local X-ray volume emissivity The model implies a higher mean black hole mass 
than is obtained from the evolution of quasars alone. 
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1987). Moreover detailed studies of the QSO X-ray lumi- 
nosity function (Boyle et al. 1994) and the source number 
count distribution have shown that QSOs are unlikely to 
form more than 50 per cent of the XRB, even at 1 keV. 

Models based on the unified Seyfert scheme, in which 
a large fraction of the emission from AGN is absorbed by 
obscuring matter (Madau et al. 1994; Celotti et al. 1995; Co- 
mastri et al. 1995) seem to work. In these models though, it is 
difficult to obtain a smooth spectrum in the 1-10 keV band 
(Matt and Fabian 1993). Some spectral features are ex- 
pected due to the matter in the vicinity producing an iron 
edge and emission line in the intrinsic spectrum. The re- 
quirement of a smooth spectrum for the XRB has been 
recently enphasized by ASCA observations which show no 
such spectral features (Gendreau et al. 1995a, 1995b), imply- 
ing that a new faint population with a very hard, smooth 
X-ray spectrum is required. 

A promising new population of sources with harder 
mean X-ray spectra discovered in deep ROSAT images 
(Hasinger et al. 1993; Vikhlinin et al. 1995; McHardy et 
al. 1995; Boyle et al. 1995; Almaini et al. 1996) is emerging 
as a significant possibility for the missing hard component 
of the XRB. At present this population can account for only 
about 10 per cent of the XRB, but the source number counts 
are still climbing at the lowest detected fluxes. Recent deep 
ROSAT studies are beginning to resolve some of the popu- 
lation into narrow emission-line galaxies (Boyle et al. 1995; 
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Griffiths et al. 1996) with remarkably high X-ray luminosi- 
ties, typically two orders of magnitude more than that of 
galaxies observed locally (Fabbiano 1989) despite compara- 
ble optical luminosities. Although it now seems clear that 
faint galaxies are emerging as a significant new X-ray pop- 
ulation, questions on the origin and nature of their activity 
still remain unsolved. 

The spectrum of the XRB, with its ~ 30keV rollover, 
requires some spectral uniformity in its constituent sources. 
Since it resembles thermal bremsstrahlung in the 3-60 
keV range, a mechanism which produces such radiation in 
galaxy nuclei would be particularly appealing. A way of 
standardising the temperature is then required. A possibil- 
ity has emerged in relation to recent discussion of energy ad- 
vection solutions (Shapiro, Lightman & Eardley 1976; Begel- 
man 1978; Rees et al. 1982; Abramowicz et al. 1995; Narayan 
1996; Chakrabarti 1996) for accretion disks by Narayan & Yi 
(1995a,b). Beginning with the work of Shapiro, Lightman & 
Eardley (1976) and Rees et al. (1982), investigations of black 
hole accretion disks at low M(< O.OlMEdd), have focussed 
on a class of optically thin solutions where the gas is sig- 
nificantly hotter than in the local Shakura-Sunyaev (1973) 
thin disk solution. Because of the poor radiative efficiency 
of the accreting gas, in the advection-dominated solution 
most of the accretion energy is stored within the gas and 
advected radially inward. The accreting plasma in this so- 
lution is two-temperature; since the ions are much hotter 
than the electrons they maintain a thick torus (which is sup- 
ported by the ion pressure). This requires the gas density to 
be sufficiently low that ion-electron coupling via Coulomb 
collisions becomes weak and cooling via synchrotron and 
bremsstrahlung radiation is not very important. In addition, 
the gas is optically thin for Compton cooling to be modest. 
For these reasons, the ion-supported torus exists only at low 
mass accretion rates (Rees et al. 1982). 

A major attraction of this class of solution is that 
the plasma attains and maintains an electron temperature 
kT c ~ 100 keV and, since the gas is optically thin, much of 
the X-ray cooling occurs through thermal bremsstrahlung. 
X-ray spectra from such sources, when integrated over red- 
shift, should plausibly resemble the XRB. 

In this paper we explore a model in which a popula- 
tion of galactic nuclei undergoes advection-dominated accre- 
tion, perhaps as their quasar activity diminishes, and so pro- 
duces the XRB spectrum. Although we centre our model on 
advection-dominated disks (ADD), any situation in which 
an optically thin, two-temperature magnetized gas occurs 
would suffice. 

In section 2 of this paper we review the physical con- 
ditions in ADD, as discussed in Narayan & Yi (1995b). We 
determine the two plasma temperatures from balancing the 
heating and cooling, and deduce the process that dominates 
the X-ray emission. In section 3, we derive the spectrum of 
the XRB from our model. We then discuss our results and 
their implications in the last section. 



2 X-RAY EMISSION FROM THE 
ADVECTION-DOMINATED DISK 

Since most of the accretion energy in advection-dominated 
flows is stored within the gas and advected radially inward 



we expect the temperature to be very high and much of the 
radiation to come out in hard X-rays. 

2.1 Physical properties of advection dominated 
disks 

Narayan & Yi (1994,1995a) have investigated the general 
properties of advection-dominated flows and derived self- 
similar and height-integrated solutions for the continuity, en- 
ergy and momentum equations. We review here those equa- 
tions from their work (and those of Stepney & Guilbert 1983 
and Zdziarski 1988) which are relevant to the present dis- 
cussion. The solutions give us quantitative estimates for the 
angular velocity, radial velocities and density of the flow as 
a function of radius R, total mass M, accretion mass M, 
and other parameters, such as the ratio of specific heats 
7, the viscosity parameter a and the fraction of viscously- 
dissipated energy which is advected, / (from to 1); a frac- 
tion (1 — /) of the energy is radiated. 

In order to obtain the temperature of the flow we need 
an equation of state and the equation for the balance of 
heating and cooling. We allow the electron temperature T e , 
and the ions temperature Ti, to be different and we take the 
gas pressure to be given by 



Pa = Ppc 3 



pkTi pkT c 



(1) 



where p is the density of the gas, c 3 the isothermal sound 
speed, /ii and /x e the effective molecular weights of the ions 
and electrons respectively, f3 the ratio of gas pressure to total 
pressure where 

Ptot = Pg+Pm, 
Pg = /3ptot, 



Pm = (1 - /3)ptot = — , 



(2) 



p m is the magnetic pressure and B the magnetic field. We 
assume the gas pressure is in equipartition with magnetic 
pressure and do not include radiation pressure. By scaling 
the mass in units of the solar mass, M = mM@, the accre- 
tion rate in Eddington units, M — ihM-E&d, where MEdd = 
Lvdd/vc 2 = 1-39 x 10 18 mgs _1 , (rj is the standard accretion 
efficiency factor), and radii in units of Schwarzschild radius, 
R = rRs, Rs = 2.95 x 10 mem we have the following rela- 
tions for p,c 2 a ,T es , (Narayan 1995b) 



= <t>(f,P)r 'cms 2 , 
P~ i'ifi P)a~ 1 m~ 1 mr~' i t 2 g cm" 
r cs = n c u T R = (p(f, I3)a^ 1 rnr^ 1/2 



p/p, c m u cm 3 , 
(3) 



where <f>, ip, ip express the dependences on /3 and / according 
to the Narayan & Yi solutions (1995b). r es is the scatter- 
ing optical depth, or is the Thomson cross section, n e the 
number density of electrons and a is the standard Shakura- 
Sunyaev viscosity parameter. 



2.2 Heating and cooling 

We determine the ion and electron temperatures in the ac- 
creting plasma by taking into account the detailed balance 
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of heating and cooling. Because almost all of the energy goes 
into internal energy, the gas becomes much hotter than in 
the thin disk solution. In particular, if the heating is adi- 
abatic, then Ti oc n 2 / 3 while T c oc n 1 ^ 3 so the ions arc 
preferentially heated. The viscous dissipation of energy acts 
primarily on the ions, which transfer little of their energy to 
the electrons. We assume that the transfer from ions to elec- 
trons occurs only through Coulomb coupling. The cooling of 
the plasma is via electrons and occurs through a variety of 
channels. 



2.2.1 Coulomb heating of electrons by ions 

The Coulomb rate of transfer of energy in the case of hot 
ions heating cooler electrons is given by (Stepney & Guilbert 
1983) 




1000 



3 m c (kTi - kT e ) 
n e maTC- 

2 Wlr, 



'K 2 {i/e e )K 2 {i/e i ) 



In A 



2(0 C + 6»;) 2 + 1 



erg cm 3 s 1 



(4) 



where InA ~ 20, is the Coulomb logarithm, the K's are the 
modified Bessel functions, and e ,0i are the dimensionless 
electron and ion temperatures defined by 



fcT c 
m e c 2 



(5) 



Electrons are mainly cooled by bremsstrahlung, cyclo- 
synchrotron and Compton cooling via interaction with soft 
photons. 



2.2.2 Bremsstrahlung cooling 

The total bremsstrahlung cooling rate per unit volume is the 
sum of electron-ion and electron-electron bremsstrahlung 



Qhi — QeiB + QccB • 



(6) 



According to Stepney & Guilbert (1983) and taking into ac- 
count the corrections to the numerical constants by Narayan 
& Yi (1995b) we adopt 

1.25n 2 aTcafm e c 2 Fi C B(0 c ) 



Qeii 



erg cm 3 s 1 



(7) 



ctf is the fine structure constant and the function F c i(8 e ) for 
both e < 1 and e > 1 is given in Narayan & Yi (1995b). 
For electron-electron bremsstrahlung we have the following 
expressions (Gould 1980; Narayan & Yi 1995b) 



< 1 



Qe" 



= n 2 c r 2 m c z a { (44 - 37r 2 )0 3/2 

x (1 + 1.10 e + el - 1.256>e /2 ) erg cm" 3 s" 1 , (8) 

0c > 1 

QocB = n 2 r 2 m o c 3 a f 240 e (ln2r/0 e + 1.28) 

ergcm~ 3 s _1 , (9) 

where r c = e 2 /m c c 2 is the classical electron radius and r\ = 
exp(— 7b) and 7e is Euler's constant ~ 0.5772. 



Figure 1. The plot shows the variations of ion temperature T; 
and electron temperature T e for three different values of m be- 
low the critical accretion rate (m cr i t < 10 -1B —see also Narayan 
&Yi 1995b). The three different m (1CT 17 -solid line, 10" 2 2 
—dotted line, 10~ 2 7 —dashed line) models give very similar tem- 
perature solutions, in particular the temperature solutions be- 
come virtually identical as rh decreases from m = 10~ 2,2 to 
10 — 2 7 . The accreting plasma is single temperature and virial at 
large radii r > 10 3 . Ti remains nearly virial, but T e saturates 
because of poor coupling between the ions and the electrons and 
due to a variety of efficient cooling mechanisms for the electrons. 



2.2.3 Cyclo- synchrotron cooling 

Due to the assumption of an equipartition magnetic field 
in the plasma, synchrotron emission becomes an important 
cooling process. The cyclo-synchroton harmonics are self ab- 
sorbed up to an energy, x t , at which self-absorption ceases 
to establish a local Planck spectrum, (Zdziarski 1988) 

C 



Xt = 6x c In 



ln- 



c 



where x t = hv t /rn e c 2 , 



27r 1/2 r es cxp(l/0 c 



30 c ' a f x c 



1/3 



(10) 



(11) 



and x c is the dimensionless cyclo-synchroton frequency for 
magnetic field in equipartition. 

To estimate the cooling per unit volume resulting from 
the cyclo-synchroton emission we adopt the expression de- 
rived by Narayan & Yi (1995). This expression has been 
obtained by equating the cyclo-synchroton emission to the 
Ray leigh- Jeans blackbody emission and by assuming that at 
each frequency v the observer sees a blackbody source with 
a radius determined by the condition v = Ut(R) so that 



_ 2tt 

Vsynch ~ ~R 



erg cm 3 s 1 



(12) 



2.2.4 Comptonization of cyclo-synchroton radiation 

Comptonization of soft photons becomes an important cool- 
ing mechanism in the inner regions of disk. Comptonization 
can be parametrized by means of the energy enhancement 
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Figure 2. The electron temperature solution is virtually independent of the model parameters. The upper curves in each plot represent 
the heating function. The lower curves the cooling function. In each of the three plots the intersection point of the respective curves 
determines the temperature solution for different values of the parameters. In (a) we vary the black-hole mass m, (b) the viscosity 
coefficient a, and (c) the fractional value of /, the advection parameter. 



factor rj (Dermer et al. 1991, Narayan & Yi 1995b). 



V = 1 + 



P(A-l) 
1 - PA 



1 - 



InP 



1 + J?i - r)2 



(*)'■ 



(13) 



where 

x — 
P = 
A = 



hv/m c c , 
1 — exp(- 
l + 49 c + 168 2 c . 



s), 



(14) 



The factor P is the probability that an escaping photon is 
scattered, while A is the mean amplification factor in the en- 
ergy of a scattered photon when the the scattering electrons 
have a Maxwellian velocity distribution. 

Most of the cyclo-synchroton radiation is emitted near 
the self-absorption energy x t given by equation (10). The 
cooling rate due to comptonization of this radiation is given 
by 



Q 



Comp 



Q 



synch 



m - m 



(15) 




log[i/(Hz)] 

Figure 3. Emission produced by accretion onto a 1O 7 M0 black- 
hole. Here we show the two major peaks in the emission corre- 
sponding to the two most efficient cooling mechanisms for the 
electrons, i.e. synchrotron and bremsstrahlung. The Compton 
emission is roughly indicated by the line between the two peaks. 



2.3 The plasma temperatures: thermal balance of 
electrons 

The net volume cooling rate for the accreting gas equals the 
sum of the three cooling rates given above, i.e. 

Q~ = Qsynch + Qcomp + Qbr- (!6) 

The energy balance of the electrons requires the net heating 
rate to be equal to the net cooling rate, this gives 

Q + = Q~. (17) 

This relation, together with equation (1) is solved to 
obtain Ti and T e , for a given m, m, /, a, and (3. In order 
to characterize the solution for advection-dominated AGN, 
we choose m = 10 7 ,a = 0.3,/ = 0.9 for and (3 — 0.5. A 
plot of the temperatures T\ and T c as a function of radius 
and for three different values of rh is given in Fig. 1. The 
ions achieve nearly virial temperature, 7] = 10 r _1 K, and 
the electrons are cooler, T e w 10 9 K for any given m below 
the critical value, rh < 0.1a 2 (Rees et al. 1982; Narayan & 



Yi 1995) for which advection-dominated solution would not 
exist, and for r < 10 3 . 

The stability of the solution for the electron tempera- 
ture T c « 10 9 K for r < 10 3 is a particularly strong prop- 
erty of the model. In Fig. 2 we plot the total cooling rate 
Q~ and the total heating rate Q + as function of temperature 
and we let the different parameters m,m,f,a, and (5 vary 
in the respective ranges of interest. These plots show that 
the temperature solution is essentially independent of model 
parameters as the energy balance solution, represented by 
the intersection point of the heating and cooling functions, 
always occurs at 6 e = 0.20 ±0.03 in the different plots. This 
is related to the fact that the slope of the cooling function 
is steeply increased by the Q~ yncb contribution just above 
6 a = 0.2 (Q~ h is the strongest function of C ; Eqns. (10)- 
(12))- 

We calculate the spectrum for the hot-ion sources 
including cyclo-synchrotron emission and bremsstrahlung 
(Fig. 3) given the model parameters described above. We 
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Figure 4. The X-ray background from sources with two- 
temperature torus disks at their centres (solid line) compared to 
Gruber's (1992) empirical formula (Eqn. 19) from 3 to 60 keV. 
We take into account the contribution of about 20 per cent from a 
power-law component of energy index ~ 0.7 mostly due to Seyfert 
galaxies. The temperature of the electrons C = 0.20, the evolu- 
tion parameter p = 2.6 and z max = 2.0. This is the best fit model. 



do not include the detailed spectrum due to Compton scat- 
tering of synchrotron photons as this is complex and not 
relevant for the present discussion. The bulk of the emission 
occurs in the high radio band close to the self-absorption fre- 
quency f t ~ 6 x 10 11 Hz. The position of the bremsstrahlung 
peak is determined primarily by T c . At the electron tempera- 
ture 8 C = 0.2, bremsstrahlung radiation dominates the hard 
X-ray emission. 



3 THE X-RAY BACKGROUND MODEL 

Applying the results reached above we now discuss the emis- 
sion spectrum of hot tori as a model for the XRB. The 
computed value of T c ~ 10 9 K is nearly independent of 
the radius r, m, rh, a, f throughout the two-temperature do- 
main, implying that such tori would have the same intrin- 
sic spectrum and that thermal bremsstrahlung would domi- 
nate the hard X-ray emission. We have calculated the spec- 
tral emissivity using the expressions for the electron-electron 
bremsstrahlung and electron-ion bremsstrahlung computed 
by Stepney and Guilbert (1983). For the evaluation of the 
electron-ion bremsstrahlung, the cross section in the Born 
approximation for a relativistic electron in a Coulomb field is 
given by Gould (1980). For electron-electron bremsstrahlung 
process, the cross-section is much more complicated and we 
have used the numerically integrated expression given in 
Stepney and Guilbert (1983). 

In order to obtain a model for the XRB, we consider the 
contribution af all sources distributed over a certain redshift 
range. These X-ray sources are not resolved as individual 
sources but they contribute to the XRB intensity. Here we 
describe the XRB as a superposition of sources from the 
local universe (z — 0) to the early universe based on the 
framework of the Robertson- Walker metric. We therefore 
expect the mean luminosity, the number density and the 



spectra of these sources to change with cosmic time. The 
total sum of X-rays from the whole population contributes 
to the XRB we observe locally. The intensity received from 
objects with local spectral emissivity jo (comoving j(z) = 
jo(l + z) p ) distributed up to redshift z max is therefore given 
by 



Ib(E)(< z u 



AtyHo 



f 

Jo 



(18) 



where go is the deceleration parameter, Ho the Hubble con- 
stant (we use go = 0.5 and Ho = 50kms~ ) and job is 
the total bremsstrahlung emissivity at the present epoch. 
Ib represents our computed XRB intensity, which we com- 
pare with the empirical fit to the XRB spectrum from 3 to 
60keV obtained by Gruber (1992) 

Ixhb = 7.877£T°- 29 exp(-.E , /41.13keV) 

keVs- 1 sr- 1 cm- 2 keV" 1 . (19) 

We plot Gruber's formula and our model of the XRB 
in Fig. 4. The fit is obtained by adding a power law compo- 
nent of energy index ~ 0.7 to our model. This takes into 
account a fraction (~ 20 per cent) contributed by AGN 
(mostly Seyfert galaxies) to the XRB in this energy band 
(the addition of a such component is not essential and good 
fits can be obtained by considering the model sources alone). 

We compare the model to the fit of the XRB by Gruber 
(1992) in terms of the fractional errors e = Ib — /xrb 
//xrb after normalizing at different photon energies. Nu- 
merically we find that the upper limit of the integral, z max ~ 
2, and p « 2.6 give the best fit to the XRB. The spectrum 
is not very sensitive to the value of the evolution parame- 
ter p, leading to residuals that are always < 1 per cent for 
2.5 < p < 3.0 in the 3-60 keV range. For p = the residuals 
increase suggesting that some evolution is required for such 
sources. 

We further examine these results by plotting contour 
levels for the maximum fractional errors associated with our 
fit for a given normalization (which directly corresponds to 
a given model volume emissivity, jom) and a given p (Fig. 5). 
We assume that acceptable fits are defined by the condition 
max(e) < 5 per cent; this condition constrains the value of 
the evolution parameter to be between 2.5 - 2.7. The model 
volume emissivity of the source population is then described 
by 

n2.6±0.1 



j(z)=j (l + zY 



(20) 



We find that the values of z max are plausibly anticorrelated 
with p, low z max implying high values for spectral evolution 
parameters for accectable fits. Relatively high evolution (due 
to both z max « 2 and p) seems to fit the peak of the XRB 
quite well as it best constrains the model parameters. The 
precise AGN contribution is not important provided it is less 
than ~ 30 per cent at 2 keV. 

3.1 The volume emissivity 

The source population which reproduces the XRB must also 
be consistent with the local values of X-ray volume emissiv- 
ity. 
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Figure 5. Contour levels for the maximum fractional error, 
max(e), calculated for each fit to the XRB obtained by varying p, 
the evolution parameter and jom is the local emissivity. jo is the 
average local value of volume emissivity determined by Carrera 
et al. (1995). The plot constrains the value of the evolution pa- 
rameter and of total contribution to the local volume emissivity, 
i.e. the normalization, to a very well defined region of parameter 
space; implying strong evolution of the sources (2.5 < p < 2.7) 
and a contribution of 10 per cent to the local volume emissivity. 

In recent years it has become apparent that X-ray emis- 
sion from galaxies can make significant contribution to the 
XRB. Direct studies of galaxy contribution using two-point 
cross-correlation between the hard (2 — lOkeV) XRB and 
various optical and infrared galaxy catalogues have been 
carried out (Lahav et al. 1993; Miyajii et al. 1994; Carrera et 
al. 1995) and have provided values for the local volume emis- 
sivity of the order of jo ~ (3.5-5.7) x 10 38 h 50 ergs -1 Mpc~ 3 
(hr,o = Ho/50 kms -1 Mpc -1 ). Extrapolation of this result 
back to z ~ 1 — 4 has lead to the conclusion that < 10 — 30 
per cent of the XRB could be produced by a population of 
galaxies if evolution is not included; some moderate positive 
evolution would lead to higher values. 

We can estimate a value for the local X-ray volume 
emissivity by integrating the model spectral emissivity over 
photon energy. In this way we obtain jom =< noLo >~ 0.6 x 
10 38 /i 5 o erg s _1 Mpc~ 3 in the 2-10 keV range implying that 
our model population would be contributing a fraction of 
0.10 — 0.15 to the local volume emissivity jo (Fig. 5). If we 
consider that field galaxies have a spatial density < n >~ 
lO^Mpc" 3 , we can determine the mean X-luminosity of 
the nearby galaxies cross-correlated with the XRB to be of 
the order of < L >^ 6x 10 40 erg s _1 (in the 2-10 keV band), 
about an order of magnitude higher than that of a typical 
spiral at the present epoch (Fabbiano 1989). Much of the 
X-ray luminosity of local spirals is however due to X-ray 
binaries. This suggests that luminous ion-supported tori are 
rare at the current epoch. 

Other determinations of volume emissivity have come 
from the studies of cross-correlation of faint galaxies cata- 
logues with unidentified deep ROSAT X-ray sources (Roche 
et al. 1995; Almaini 1996; Treyer & Lahav 1996) as it 
has emerged recently that faint galaxies could be impor- 
tant contributors to the XRB (see Introduction). The es- 
timated value of the volume emissivity for a galaxy sam- 



ple at a median redshift of z = 0.27, is jo ~ (2.3) x 
10 38 /i5 ergs _1 Mpc~ 3 in the 0.5 to 2keV range (Almaini 
1996, Treyer & Lahav 1996) . The value we deduce from 
our model for this energy band is jom =< noLo >~ 1.2 x 
50 erg s Mpc where we have considered a sample of 
model sources at the same redshift. Almaini (1996) has con- 
strained the evolution of the X-ray emissivity with redshift, 
obtaining evidence that X-ray emission from such galaxies 
evolves as strongly as AGN. These results are in good agree- 
ment with our model (p ~ 2.6). 

A property of these X-ray emitting galaxies is that the 
X-ray luminosity, L x ^ 10 42 ergs -1 (Boyle et al. 1995; Grif- 
fiths et al. 1995; Almaini 1996), is typically two order of mag- 
nitude more than that of galaxies observed locally. If our 
model sources had this mean luminosity their number den- 
sity would be of the order of < n >~ 10~ 4 - 10~ 5 Mpc~ 3 . 
These values of < no > constrain the spatial density of the 
model population to be very close to that of other AGN 
at ~ 1/100 the spatial density of local field galaxies. This 
determination is strengthened by the fact that the sample 
analyzed by Roche et al. (1995) and Almaini et al. (1996) 
goes deeper in redshift and possibly decreases the errors as- 
sociated with the uncertain evolutionary properties of the 
model galaxies. Both the lower luminosity objects observed 
locally and the higher luminosity population at z = 0.3 are 
in agreement with our ion-torus scenario in which evolution 
is important, the accretion rate is typically low, M < 1 
per cent of JVfEdd and the black-hole mass is of the order of 
M — 10 7 — 10 9 M Q (the hard-spectrum is essentially inde- 
pendent of black-hole mass) . 

3.2 Radio emission 

The bulk of the emission from the ion-supported tori, occurs 
at radio and submillimeter band near the self-absorption 
frequency v t ~ 6 x 10 1 Hz (see Fig. 3). Directly from Fig. 3 
we determine a mean radio flux density for the ion-supported 
tori model at ~ 10 9 Hz and compare it with the results 
of deep VLA images of radio sources. Hamilton & Helfand 
(1993) have found some corrrclation of faint radio sources 
at 20 cm with XRB fluctuations and associated them with 
the unresolved XRB. We find that our model predicts a flux 
density of the order of « 0.3 mjy at 10 9 Hz (near 20 cm) for 
a source at z = 0.3. The derived value of the flux density is 
consistent with the lower end of the range of flux densities 
found by Hamilton & Helfand (1993). 

Radio surveys at high frequencies (closer to the peak of 
the self-absorption frequency) are potentially best for select- 
ing ion-supported tori. They could confirm our model if the 
radio sources correlate with faint, non-AGN, X-ray sources 
and have rising radio spectra. 



4 SUMMARY AND DISCUSSION 

The diffuse XRB in the few-lOOkeV band can be well ex- 
plained as due to a population of apparently normal galaxies 
undergoing two-temperature advection-dominated accretion 
via an ion-supported torus. 

Here, we have considered a physical model for the de- 
scription of the high energy emission and the cosmologi- 
cal evolution of such sources. We find that a population of 
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sources with X-ray emission due to an optically-thin thermal 
plasma emitting bremsstrahlung at a temperature 8 a ~ 0.2 
satisfies the background constraints. We show that the tem- 
perature requirement is plausibly achieved in ion-supported 
tori. 

The spectral emissivity of such sources is most likely 
to be associated with (some) spectral evolution and the ac- 
ceptable evolution parameters are roughly consistent with 
those of AGN. The predicted population of sources also 
satisfies the constraints of the local volume emissivity. In 
particular we have found that the predicted emissivity con- 
tributes a small fraction to the local volume emissivity in 
the 2-10 keV band. For the deeper ROSAT sample, the spec- 
tral emissivity extrapolated in the 0.5-2 keV more closely re- 
sembles the value obtained with our model, suggesting that 
advection-dominated disks might be related to the activity 
of X-ray luminous, narrow-line galaxies. 

We note that the low radiative efficiency in the 
advection-dominated scenario may indicate that most galax- 
ies might harbour a massive black hole as a remnant of 
earlier more active quasar epoch (see also Fabian & Rees 
1995). In particular, if we suppose that QSOs contribute 
about 20 per cent of the XRB at 1 keV and have power-law 
spectrum of energy index ~ 1, we can calculate a ratio of 
the total emission due to QSOs to that of ADD sources of 
Eqso/Eadd ~ 1/4. If we now assume that the radiative 
efficiency r\ of QSOs and ADD tori to be about 10 per cent 
and 0.1 per cent respectively, we derive an upper limit for 
the ratio of the mass accretion rates, i.e., 

Mqsq < 1 / ??add \ < 1 
Madd ~ 4 I 0.1 ) ~ 400' 

This leads directly to the conclusion that the mean mass of 
'dead' black holes might be up to a factor of 100 greater 
than that estimated previously from measurements of the 
luminosities. Moreover, given that most of the XRB could 
be due to a population of ion-supported tori, we can find 
limits for the total mass density, p, of black-holes distributed 
in local galaxies (Soltan 1982 or Fabian & Canizares 1988). 
The energy density due to accretion is given by e = prjc 2 
which is equivalent to the the total energy emitted by all of 
the model sources in a unit comoving volume; i.e. the model 
XRB emission. If we take the spatial density of host galaxies 
to be « 10 -3 Mpc -3 , we find that most normal galaxies 
should contain a central black hole with a mass w 10 7 M Q . 
Lower space densities lead of course to higher mean masses. 
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